Trapping of light pulses in ensembles of stationary A atoms 
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We present a detailed theoretical description of the generation of stationary light pulses by stand- 
ing wave electromagnetically induced transparency in media comprised of stationary atoms. We 
show that, contrary to thermal gas media, the achievable storage times are limited only by the 
ground state dephasing rate of the atoms, making such media ideally suited for nonlinear optical 
interactions between stored pulses. Furthermore, we find significant quantitative and qualitative dif- 
ferences between the two types of media, which are important for quantum information processing 
schemes involving stationary light pulses. 
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I. INTRODUCTION 



The coherent transfer of quantum states between light 
and atoms has been the subject of much research, both 
experimentally and theoretically, motivated by potential 
applications in quantum computing, quantum cryptog- 
raphy and teleportation. While the transfer of quan- 
tum states from light to a single atom can in principle 
be achieved by cavity QED techniques [![, the required 
strong-coupling regime is experimentally very difficult to 
reach. To overcome this difficulty the use of atomic en- 
sembles, rather than single atoms, has been proposed 
@, H, HJlland implemented for storage of classical light 
pulses [3 3 j an d recently storage of non-classical pulses 
has been demonstrated [1, One such light storage 
scheme is based on electromagnetically induced trans- 
parency (EIT) [To[| in ensembles of A-type atoms Q. 
While the storage and retrieval of light pulses with this 
scheme has been demonstrated utilizing many different 
types of atomic ensembles, including Bose-Einstein con- 
densat es f§j ar id thermal gasses Q as well as solid state 
media [111 Il2| , the non-trivial manipulation of, and in- 
teraction between, stored pulses is hampered by the in- 
herent trade-off between storage time and field ampli- 
tude. A step towards overcoming this problem was taken 
with the suggestion of using standing wave fields to cre- 
ate a periodic modulation of either the dispersive [l3| or 
the absorptive [lij properties of the medium, inducing 
a photonic bandgap [T5| and creating a stationary light 
pulse. Schemes to implement a controlled phase gate us- 
ing these techniques have been proposed using either the 
dispersive or the absorptive [l?} grating technique, 
but both are still hampered by a trade-off between stor- 
age time and field amplitude. For the latter case, only 
thermal gas media have been considered, and a detailed 
theoretical treatment of this case is given in [18]. This 
theory, however, does not apply to media comprised of 
stationary atoms such as ultra cold gasses or solid state 
media [19]. In this article we present a detailed theoret- 
ical treatment of the creation of stationary light pulses 
by the absorptive grating technique for media comprised 
of stationary atoms. We find that the loss of excitations 



inherent to the thermal gas case is absent for station- 
ary atoms, making such media ideally suited for the kind 
of non-linear optical interactions envisaged in [17]. Fur- 
thermore we find interesting quantitative and qualitative 
differences between the thermal gas and ultra cold gas 
cases when quasi-standing wave coupling fields are con- 
sidered. These differences are imp ortant for the proposed 
controlled phase gate scheme [17| in stationary atom me- 
dia. 

In Sec.|n]we present a detailed account of our theory of 
stationary light pulses in media comprised of stationary 
atoms and compare the results to the thermal gas case. 
The theory is complemented in Sec. IIIII by a calculation 
of non-adiabatic corrections. A summary of our results is 
provided in Sec. II VI Appendix [Al contains a brief review 
of the theory of stationary light pulses in thermal gas 
media [l8|], reformulated in terms of polariton fields, used 
for comparison with the stationary atom case. 



II. STANDING WAVE POLARITONS IN 
ENSEMBLES OF STATIONARY ATOMS 



We consider an ensemble of N non-moving A atoms 
interacting with probe and coupling lasers propagating 
parallel to the z axis. The two lower states \b) and |c) 
of the atoms (see Fig. ) are assumed to be nearly de- 
generate, such that the magnitude of the wave vectors of 
the probe and coupling lasers can be considered identical 
{hp — k c — /c) . 

The Hamiltonian for the N atom problem is 
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H J 



HI 



H v 



(1) 



where Hp and Ha describe the free electromagnetic field 
and the atoms, Hi describes the interaction of the atoms 
with the probe and coupling fields, and Hy describes the 
interaction with the vacuum field modes. The individual 
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FIG. 1: The 3-level A atom. The quantized probe field couples 
the ground state |6) to the excited state \ a), while the classical 
coupling field couples the metastable state |c) to |o). 



terms are given by 
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We introduce slowly varying field operators for the elec- 
tromagnetic field. Since we are allowing for standing 
wave fields, we write the field operator as a superposi- 
tion of two traveling wave fields propagating in opposite 
directions 
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in the rotating wave approximation are 
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where 7 — j b + "f c is the decay rate of the excited state 
I a) into the two lower states. The complex decay rates 
r M „ are given by 

(7) 
(8) 
(9) 

where ip V are the dephasing rates of the respective co- 
herences, S Pt c are the one-photon detunings of the probe 
and coupling lasers, respectively, and A is the two-photon 
detuning. We have also assumed that the coupling field 
can be treated as a classical field with Rabi frequency Sl c 
given by 
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(10) 



In the following we shall disregard the noise operators 
F^y since we will be considering the adiabatic limit. 
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where the operators E pc are given by 

Ep tC (z, t) = E+ C (z, t)e ikz + E- C (z, t)e- lkz 



(3) 



(4) 



The field operators Ep C are the slowly varying field oper- 
ators for the forward and backward propagating compo- 
nents of the probe and coupling fields with carrier fre- 
quencies Wj,,c> and e p c are the respective polarization 
vectors. 

We define continuum atomic operators <7 Mt , by sum- 
ming over the individual atoms in a small volume V, and 
introduce slowly varying atomic operators a^ v defined by 
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A. Weak probe approximation 

In order to solve the propagation problem, we assume 
that the probe field is weak compared to the coupling field 
and that the probe photon density is small compared to 
the atomic density. In this case the Heisenberg-Langevin 
equations can be solved perturbatively. To first order 
in the probe field amplitude, the relevant Heisenberg- 
Langevin equations are 
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Combining equations (jlip we can obtain a differential 
equation for a bc 
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B. Adiabatic limit 



Notice that the operators defined by ((5]) are slowly vary- 
ing in time, but not in space. 

The Heisenberg-Langevin equations for these operators 



In order to solve equation (|T2|) , we consider the adia- 
batic limit in which the fields vary slowly in time. Intro- 
ducing a characteristic timescale T of the slowly varying 
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operators, we expand Obe in powers of (7b a T) . To ze- 
roth order we find 



(13) 



Inserting this expression into (|lla[) . we find an expression 
for a^a valid in the adiabatic limit 



1 / d \ (g P E p 

a ba = -T77- 1 be + — 



dt \ n. 



(14) 



By inserting the field decomposition Q into the adia- 
batic expression for <Jba (|14|) we get 



-g P (r bc + j|) 



+ 2|«+||k-| cos(2fcz + 4>)) 
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where we have also introduced the time dependent total 



Rabi frequency f2(t) = y |fi^| 2 + \Q C | 2 and the ratios 

= -fi- which are assumed to be constant. The phase 
angle <j> is defined by the relation 
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Polariton field 



(16) 



We now introduce a dark-state polariton (DSP) field 
analogous to the DSP field defined in Q 

E±(z,t) = coB0(t)¥ ± («,t), (17) 

where the angle 8 is given by the total coupling laser Rabi 
frequency through 



tan0(i) 



fi(t) 



(18) 



By inserting the definition (p~7|) of the DSP field into fT5 
we get 



be 



+ 2\k+\\k~\ cos(2fcz + 4>)) mqs 
x (smd(^ + e lkz + y-f- lkz ) 



N r <Tba 



To derive wave equations for the components of the DSP 
field, we need to expand the optical coherence aba in 
spatial Fourier components. We do this by inserting the 
Fourier series 



1 _ ao 
1 + y cos x 2 



= — + y^a n cos(rtx), 



(20) 



where y = 2|k + ||k~| and x — 2kz + <fi, into (jT9j) . Note 
that y < 1 which guarantees the existence of the Fourier 
series except in the case of a standing wave coupling field 
(y = 1). Fortunately, we can treat this case successfully 
by considering the limit y — ► 1 at the end of our calcula- 
tion. 

Inserting the Fourier series into (119p we find 
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From (f2Tj) we see that <Jb a can we written as 
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(22) 



To derive a set of wave equations for the polariton field 
components, we need to calculate the components a^ 1 
and cr^ 1 of the expansion (|22|) which we label cr^j for 
brevity. These components are given by 

= ^ (r bc + [sin% * + + axe***")] 



(23a) 

+ )]• 

(23b) 



We see that we only need to calculate the first two Fourier 
coefficients ao and a\ of the expansion (|2ip . These are 
given by 
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(24a) 
(24b) 



Inserting the adiabatic expression (|23|) for cr^ into the 
wave equations for the probe field components 



(25a) 



we obtain a set of coupled wave equations for the DSP 
field components 
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where we have introduced a new angle 9' defined by 



tan0' 

as well as the constant 

s = — 

a 



tan 9, 



(27) 



(28) 



Since we are considering the adiabatic limit in which the 
coupling field Rabi frequency changes slowly in time, we 
shall neglect the last term on the rhs. of (|26p in the fol- 
lowing. 



D. Low group velocity limit 

In the experimentally relevant low group velocity limit 
cos 2 9 <C 1, the wave equations (f2"6")) take the simpler form 
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in the case where |k + | > \k |. In the opposite case, 
|« + | < \k~\, the wave equations are 



be 



dt 
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We have introduced the group velocity v g — c cos 2 9 
in equations (I29|) and (|30|) . and have also made use of 
the fact that in the low group velocity limit, cos 2 9' ~ 
^1 — y 2 cos 2 0. 



E. Initial conditions 

We shall consider the same kind of experiment as in 
in which a probe pulse, propagating under the in- 
fluence of a copropagating traveling wave coupling field, 
is stored in the medium and subsequently retrieved by a 
standing wave coupling field with > |k — |. 

Assuming that the standing wave coupling field is 
switched on at t = 0, we need to find the initial condi- 
tions for the two components of the DSP field 4 ,± (z,0). 
The initial condition for the Raman coherence is 



(31) 



where ^(z,0) is a known function of z determined by 
the DSP field prior to switching on the standing wave 
coupling field. Using (fl~3|) and the definition of the DSP 
field in the standing wave case (|17|) , along with the initial 
condition (|3"Tj) . we get 

¥(z,0) (n + e lkz + n- e - tkz ) = ^ + (z,0)e lkz +^-(z,0) e - lkz . 

(32) 

From this expression we see that the initial conditions for 
the components of the DSP field 4 ,± (z,0) are 

*+(z,0) = k+*(z,0), y-(z,0) = K~V(z,0). (33) 

With the initial conditions (l33l). we find the solution 



* + (M) = V 







,-.+ |2 



*(z-0r(t),O) + 1 



¥(* + j8r(i),0) 



-r bc t 



- (3r(t), 0) + *(z + /3r(t), 0) 
where /3 = v / |K+p(|^+p _ | K -|2) and r ( f ) = ccos 2 0( t ')dt'. 



(34a) 
(34b) 



5 



Probe retrieval by a standing wave coupling 
field 



To determine the solution for a standing wave coupling 
field, we let — ► A= in the solution In this limit 

we find the solution 



*-(z,t) = -^=*(z,0)e- r >>° t . 



(35a) 
(35b) 



As an example, we take the initial condition for the 
DSP field to be $(z, 0) = * exp(-(z/L p ) 2 ), where L p is 
the characteristic length of the stored probe pulse. The 
polariton amplitude 'J'o is related to the initial probe field 
amplitude Eq by = Eq/ cosOq, where 8q is determined 
by the Rabi frequency of the traveling wave coupling field 
prior to storage. 

The components of the retrieved probe field found from 
([17]) are 

E+(z,t) = 1 £^M £ ; exp(-(z/L p ) 2 )e- r - t , (36a) 
V2 cos W 

E;(z,t) = J-^M Eo exp(~(z/L p ) 2 )e~ r ^. (36b) 
V2 costfo 

In Fig. [5] we compare the retrieval of an initially stored 
probe pulse by a standing wave coupling field in the ther- 
mal gas and ultra cold gas cases. The time dependence 
of the angle 9 is assumed to be given by cos 2 6(t) — 
cos 2 6*o tanh(i/T s ) for t > 0, where T s is the characteris- 
tic switching time. For simplicity, we have assumed zero 
Raman dephasing (r^c = 0) and taken the characteris- 
tic length of the stored probe pulse to be L p = v gi oT s , 
where v g ^ — ccos 2 #o- The probe field photon density 
averaged over many wavelengths \E+\ 2 + |-Ep~| 2 , in units 
of the photon density prior to storage |-Eb| 2 , is plotted as 
a function of z in units of L p and t in units of T s . In both 
the stationary atom case and the thermal gas case we see 
that the stored probe pulse is revived into a stationary 
probe field, but we note that in the stationary atom case, 
the diffusive broadening of the probe field, evident in the 
thermal gas case, is absent. 

The solution for thermal gas media is based on the the- 
ory in [l8| , which is reviewed briefly in appendix [AJ The 
medium is characterized by the absorption length in the 
absence of EIT l a = 0.1x L p , which roughly corresponds 
to the conditions in uM. 
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(a) Stationary atoms 



(b) Thermal gas 



FIG. 2: (color online) Retrieval of a stored probe pulse with a 
standing wave coupling field. The probe field energy density, 
in units of ^r-\Eo\ 2 , is plotted for a medium comprised of 
stationary atoms (a) and thermal atoms (b) as a function of 
z in units of the pulse length L p , and t in units of the switching 
time T 3 . The absorption length of the media is taken to be 
l a = 0.1 x L p . 



of the revived probe pulse in the direction of the stronger 
of the two coupling field components. 

In the ultra cold gas case considered here, we find from 
the solution that the revived probe pulse instead 
splits into two parts. A stronger part which propagates 
in the direction of the stronger of the coupling field com- 
ponents, and a weaker part which propagates in the op- 
posite direction. 

Fig. [3] shows the solution (J5H) with the same ini- 
tial conditions as in Fig. [21 but with k + = V0.55 and 
kT = V0.45. Fig. 2] compares the retrieval of a stored 
probe pulse by a quasi-standing wave coupling field in 
thermal and ultra cold gas media. The splitting of the 
revived probe pulse is clearly evident in the cold gas case, 
indicating a qualitative difference between the thermal 
gas and the ultra cold gas cases. The cause of this differ- 
ence is the coupling to the high spatial-frequency compo- 
nents of the Raman coherence Obc in the ultra cold gas 
case. This splitting of the probe pulse is very important 
when considering various schemes for interacting pulses. 
In the phase-gate proposal of Andre et al. [l7j . a small 
imbalance in the two components of the coupling field is 
used to propagate a quasi-stationary light pulse across 
a stored excitation in a thermal gas medium. As is ev- 
ident from Fig. 0] this scheme would not work in media 
comprised of stationary atoms, since a large part of the 
revived probe field would then propagate in the wrong 
direction. 



H. Calculation of the Raman coherence 



G. Probe retrieval by a quasi-standing wave 
coupling field 

We shall now study the situation in which the probe 
field is retrieved by a quasi-standing wave coupling field. 
In the previous section, we saw that in the thermal gas 
case a quasi-standing wave coupling field leads to a drift 



To calculate the Raman coherence of the atoms, we 
use the zeroth order expression (| 13[) for a\y C 



(37) 



Sir 



By inserting the decompositions of the probe and cou- 
pling fields, as well as the definition (|17p of the DSP 
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For the rapidly varying components of the Raman coher- 




FIG. 3: Retrieval of a stored probe pulse with a quasi- and 
standing wave coupling field (k + — \/0.55, K ~ = V0.45). Fig- 
ures (a) and (b) show the polariton amplitudes in units of 

*I/o as a function of z and t, in units of L p and T s , respectively. V N r <r bc = (43) 



' 0.-2, 2 ' 0.-2, 2 





(a) Cold gas 





(b) Thermal gas 



where, in both cases, n > 0. 

In the case of a perfect standing wave coupling field, 
only the dc component of the Raman coherence is present 
which is given by 



N,a£\z,t) = -fnn6(t)*(z,0)e 



-r bc t 



(44) 



FIG. 4: Retrieval of a stored probe pulse with a quasi- 
standing wave coupling field. The probe field energy density, 
in units of \ Eo j 2 , is shown for both the ultra cold gas case 
(a) and for the thermal gas case (b). Parameters are the same 
as in Fig. [3] 



field, we get 

'iW = - sin 9 I ^ : _ v _l J . (38) 



In the quasi-standing wave case, the rapidly varying com- 
ponents of the Raman coherence o~¥t with negative val- 
ues of n attain a small but non-vanishing value, becom- 
ing progressively smaller with decreasing n. The rapidly 
varying components of the Raman coherence with posi- 
tive values of n all vanish. An asymmetry in the Raman 
coherence is to be expected, since neither the coupling 
field nor the revived probe field is symmetric in z. 



Inserting the solution (134|) into this expression, we find 
by a binomial expansion 



III. NON-ADIABATIC CORRECTIONS 



r K(Jbc = -\sme\ >!'(: >Y.()i + *C + h-M) 



+ 



^[¥(*-/3r,0)-¥(* + /3r ) 0)] (39) 



n=0 



From this expression we see that the Raman coherence 
can be written as 



In [20j it was shown that the finite length of the probe 
pulse leads to a broadening of the pulse envelope due to 
dispersion. In this section we shall investigate the same 
effect in the standing wave case and show that the disper- 
sive broadening vanishes in the case of a pure standing 
wave coupling field. 

Our starting point is the differential equation (fl2|) for 



the Raman coherence <7{, c 
find 



To first order in (jbaT) we 



(40) 



n— — oo 



where the dc component is 
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*(s-/3r(t),Q) 



(3 



9(z + 0r(t),O) 



(41) 



-T ba t 



where we have assumed Tbc — to simplify the calcula- 
tions. Inserting this expression into (|lf ap and introduc- 
ing the DSP fields defined in (jTTJ) , we get 
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'N r a ba = 



sin 9 



d 



iO(l + 2|/t+||«;-|cos(2fez + ^)) dt 
T ba sin 9 tan 2 9 d 2 

9p N r 01(1 + 2\k+\\k-\ cos(2kz + (j))) 2 dt2 

I 



-\- ikz 



* e 



-ikz\ 



(46) 



where we have assumed that the coupling laser Rabi fre- 
quency changes slowly enough to set 9 = in the equa- 
tions. 

As in section|n]we need to find the Fourier components 
err - To do this we apply the Fourier series ([2D]) and we 
introduce 



The Fourier coefficients 00,1 have already been calculated 
and are given by (|24|) , while the Fourier coefficients c?o,i 
are given by 



1 



(f + ycosx) 2 



do 
2 



d n cos(nx) 



(47) 



where, as before, y = 2\k + \\k | and x — 2kz- 
ing the Fourier series into (j4"6"]) . we get 
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(49b) 



We now insert the expressions (148)) into ([25]) to obtain a 
set of coupled wave equations for the DSP fields 
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(50a) 
(50b) 



where we have introduced the constants Once again we consider the low group velocity limit 

cos 2 fl<l. With this approximation the wave equations 
s= Sl s> = S± s " = Sl, (5i) simplify to 

ao ao ao 



ccos 



, 99+ 

dz 



o ir t> 
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-se 
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~dt 



9 2 p Nr 
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(52a) 
(52b) 



To the same order of approximation, we can replace the second time derivatives of the DSP fields with the second 
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time derivative of the zeroth order solution. Differenti- 
ating both sides of (|5l2|) with respect to t, and discarding 
derivatives of order greater than two, we get 



d 2 ^+ 



dz at 
-ccos 6 - — — —se 



,d 2 ^+ 



(53a) 
(53b) 



dt 2 " dz at dt 2 ' 

where we once again assume that the coupling laser Rabi 



frequency changes slowly. Using (|52p we can solve for the 
second time derivatives of the DSP field. We find 



d 2 ^ 



i-y 2 

l-s 2 



ccos 



dz 2 



(54) 



where we exploited the fact that in the low group ve- 
locity limit cos 2 0' ~ yl — y 2 cos 2 9. Inserting (|54jl into 
(|5^|) . and assuming that |k + | > | « | , the coupled wave 
equations take the form 



d^ 
~dt 



xl +|2 5 * + + - 

+ \K \ 1> — = K K V, 



9*- \ K +\ 2 l a v g d 2 



J9 dz 



\n + \ 2 v g 



dz ^1^2 Oz 2 

\n+\ 2 l a v g d 2 



9 dz ^/T^y 2 dz 2 



(|k+| 2 *+ - k + k~*9~) 



r 



To solve the coupled wave equations (|55|) we proceed ^/ ± (2;,t) — > VE ,± (<7, t), and find the solution 
by Fourier transforming with respect to z, such that 



(55a) 
(55b) 



where 



and 



* + M)= {bV-(q,0) - (\ K +\ 2 - d)V + (q,0)) e^q\+r(t)) 



+ ({\k+\ 2 + d)V + (q,0) - 6*~(?,0)) exp(i ? A_r(i)) 
*"(?,*) = ^ f (-&**+(?, 0) + (|k+| 2 + d)*-(?,0)) exp( l(? A + r(t)) 
+ (&*^+(g,0) - (|k+| 2 - d)§-(?, 0)) exp(igA_r(t)) 

I 



(56a) 



(56b) 



K + K-*(l-iq£), A d 



\K + \ 2 lg 



(57) 



(58) 



From this solution it is clear that the broadening of the 
pulse envelope due to dispersion is absent in the case of a 
pure standing wave coupling field. The effect is present 
in the case of a quasi-standing wave coupling field. If we 
consider the limiting case of a traveling wave coupling 
field (k - — > 0), we find the same dispersion term in the 
wave equation (|55p that is given in [20j. 



vV+l 2 (\k+\ 2 - |k-| 2 ) - |k+| 2 |k-| 2 £V. (59) 



Inserting the initial conditions (|33| for the DSP field 
and considering the limit k — > -4= , corresponding to a 
pure standing wave coupling field, the solution becomes 



(60a) 
(60b) 



IV. SUMMARY 

In this article we have presented a detailed theoreti- 
cal treatment of stationary light pulses in media com- 
prised of stationary atoms, such as ultra cold gasses and 
solid state media. We found that contrary to the thermal 
gas case, the achievable trapping time is limited only by 
the Raman dephasing rate of the atoms and such media 
are thus ideally suited for the kind of nonlinear opti- 
cal interactions envisaged in fl7T |. It was also shown 
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that the behavior of the probe pulse when employing 
quasi-stationary coupling fields is significantly different 
for moving and non-moving atoms. This fact must be 
taken into account when considering schemes for inter- 
acting pulses. Although, to the best of our knowledge, 
no experiment with stationary light pulses in ultra cold 
media has yet been reported, several experiments on nor- 
mal E1T and light storage have been performed with ul- 
tra cold gasses [(| and solid state media [ll|, [HJ . These 
experiments have also demonstrated the possibility of us- 
ing beam geometries other than copropagating probe and 
coupling lasers. We therefore expect that the experimen- 
tal demonstration of stationary light pulses in such media 
is within present day capability. 



gas case presented in [18j. It is argued that the motion 
of the atoms in a thermal gas causes a rapid dephasing 
of the spatially rapidly varying components of the Ra- 
man coherence and it is therefore assumed that only the 
n = component in the expansion (|40[) is non- vanishing. 
Consequently, the only non- vanishing components of the 
optical coherence in the expansion is the n = — 1 
and n = terms. With this approximation, the relevant 
Heisenberg-Langevin equations for the slowly varying op- 
erators are 



'ba 



'ba 



ifl~a bc - r ba a^ a , 



(Ala) 
(Alb) 
(Ale) 
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APPENDIX A: STANDING WAVE POLARITONS 
IN THERMAL GASSES 

As shown in Sec.HJJthe behavior of the stationary light 
pulses depends critically on whether the EIT medium is 
comprised of stationary or moving atoms. In this section 
we present a brief review of the theory for the thermal 



As shown in [l8[ these equations can be solved approx- 
imately by adiabatically eliminating the optical coher- 
ences cr^j and making an adiabatic expansion of (|Alc[) . 
The resulting expressions for the components of the op- 
tical coherence is then inserted into the wave equa- 
tions (|25p . Contrary to [l8|, which deals directly with 
the probe field operators E^, we introduce the polari- 
ton field defined by (|17[) which enables us to treat time 
dependent coupling fields in a consistent manner. To fa- 
cilitate the solution of the resulting wave equations, sum 
and difference normal modes defined by 



(A2a) 
(A2b) 



are introduced. In the case of an optically thick medium, 
the difference mode can be adiabatically eliminated, re- 
sulting in a diffusion equation for the sum normal mode 



dt 



Ik | 2 )ccos 2 



} 99 S 

dz 



4k + | 2 |K-| 2 Z a ccos 2 



dz 2 



T bc sin 2 0* « 



(A3) 



where we have assumed zero probe field detuning (S p = The solution of (|A3|) . subject to the initial conditions 
0). The difference normal mode is given by (f33|) . is the basis for the comparison between thermal gas 

media and stationary atom media presented in Sec. [Til 

* D = -2K+K-Z a ^. (A4) 
oz 
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